Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

Wetlands and low-gradient topography are associated with longer hydrologic transit times in Precambrian Shield headwater catchments

The estimation of hydrologic transit times in a catchment provides insights into the integrated effects of water storage, mixing dynamics, and runoff generation processes. There has been limited effort to estimate transit times in southern boreal Precambrian Shield landscapes, which are characteristically heterogeneous with surface cover including till, thin soils, bedrock outcrops, and depressional wetland features that play contrasting hydrologic roles. This study presents approximately 3.5 years of precipitation and streamflow water isotope data and estimates mean transit times (MTTs) and the young water fraction (p_y) across six small catchments in the Muskoka-Haliburton region of south-central Ontario. The main objectives were to define a typical range of MTTs for headwater catchments in this region and to identify landscape variables that best explain differences in MTTs/p_y using airborne light detection and ranging and digital terrain analysis. Of the transit time distributions, the two parallel linear reservoir and gamma distributions best describe the hydrology of these catchments, particularly because of their ability to capture more extreme changes related to events such as snowmelt. The estimated MTTs, regardless of the modelling approach or distribution used, are positively associated with the percent wetland area and negatively with mean slope in the catchments. In this landscape, low-gradient features such as wetlands increase catchment scale water storage when antecedent conditions are dryer and decrease transit times when there is a moisture surplus, which plausibly explains the increases in MTTs and mean annual runoff from catchments with significant coverage of these landscape features.     

Changes in hillslope hydrology in a perched,shallow soil system due to clearcutting and residual biomass removal

Sustainable fuels legislation and volatility in energy prices have put additional pressures on the forestry sector to intensify the harvesting of biomass for “advanced biofuel” production. To better understand how residual biomass removal after harvest affects forest hydrology in relatively low slope terrain, a Before-After-Control-Impact (BACI) study was conducted in the USDA Forest Service's Marcell Experimental Forest, Minnesota, USA. Hydrological measurements were made from 2010–2013 on a forested hillslope that was divided into three treatment blocks, where one block was harvested and residual biomass removed (Biomass Removed), the second was harvested and residual biomass left (Biomass Left), and the last block was left as an Unharvested Control. The pre-harvest period (2 years) was 2010–11 and post-harvest (2 years) was 2012–13. Water table elevation at the upslope and downslope position, subsurface runoff, and soil moisture were measured between May–November. Mixed effect statistical models were used to compare both the before-after and “control” treatment ratios (ratios between harvested hillslopes and the Unharvested Control hillslope). Subsurface runoff significantly increased (p < .05) at both harvested hillslopes but to a greater degree on the Biomass Left hillslope. Greater subsurface runoff volumes at both harvested hillslopes were driven by substantial increases during fall, with additional significant increases during summer on the Biomass Left hillslope. The hydrological connectivity, inferred from event runoff ratios, increased due to harvesting at both hillslopes but only significantly on the Biomass Left hillslope. The winter harvest minimized soil disturbance, resulting in no change to the effective hydraulic conductivity distribution with depth. Thus, the observed hydrological changes were driven by increased effective precipitation and decreased evapotranspiration, increasing the duration that both harvested hillslopes were hydrologically active. The harvesting of residual biomass appears to lessen hydrological connectivity relative to leaving residual biomass on the hillslope, potentially decreasing downstream hydrological impacts of similar forestry operations.     

Symplectic Mappings of Co-orbital Motion in the Restricted Problem of Three Bodies

The dynamics of co-orbital motion in the restricted three-body problem are investigated by symplectic mappings. Analytical and semi-numerical mappings have been developed and studied in detail. The mappings have been tested by numerical integration of the equations of motion. These mappings have been proved to be useful for a quick determination of the phase space structure reflecting the main characteristics of the dynamics of the co-orbital problem.     

Optical properties of solid and liquid sulfur at visible and infrared wavelengths

Knowledge of the optical constants of elemental sulfur has potential applications to Venus, Jupiter, Io, Amalthea, and the Earth. The real part, n, of the index of refraction of liquid sulfur (at 133°C) and of solid orthorhombic sulfur (at 25°C) for the wavelength range 0.4–2.0 μm were measured ellipsometrically. The imaginary part, k, of the refractive index of liquid sulfur was obtained by transmittance measurements at the same temperature and wavelength range. The reflectance of semi-infinite slabs of solid and liquid sulfur is calculated using the measured n and k values. We confirm that sulfur melts on Io would be classified as “black” by the Voyager imaging system.     

A focal plane detector design for a wide-band Laue-lens telescope

The energy range above 60 keV is important for the study of many open problems in high energy astrophysics such as the role of Inverse Compton with respect to synchrotron or thermal processes in GRBs, non thermal mechanisms in SNR, the study of the high energy cut-offs in AGN spectra, and the detection of nuclear and annihilation lines. Recently the development of high energy Laue lenses with broad energy bandpasses from 60 to 600keV have been proposed for a Hard X ray focusing Telescope (HAXTEL) in order to study the X-ray continuum of celestial sources. The required focal plane detector should have high detection efficiency over the entire operative range, a spatial resolution of about 1mm, an energy resolution of a few keV at 500keV and a sensitivity to linear polarization. We describe a possible configuration of the focal plane detector based on several CdTe/CZT pixelated layers stacked together to achieve the required detection efficiency at high energy. Each layer can operate both as a separate position sensitive detector and polarimeter or work with other layers to increase the overall photopeak efficiency. Each layer has a hexagonal shape in order to minimize the detector surface required to cover the lens field of view. The pixels would have the same geometry so as to provide the best coupling with the lens point spread function and to increase the symmetry for polarimetric studies.     

Problems in data compilation

Concepts involved in the estimation of target quantities and other adjustment parameters are critically discussd. We point out that one can find more accurate precepts for the reduction of data by utilizing all available constraints on all available data in the derivation of the reduction precepts. We introduce a measure for theefficiency of a set of adjustment parameters such that adjustments carried out using different precepts can be objectively compared. Finally, having applied our suggestions to a specific problem, we show that we have obtained estimates of a set of target quantities (in our case, star positions and proper motions) which have smaller formal errors than estimates of the same target quantities derived from the same input material but following traditional procedures.     

Photometric lightcurves of 21 asteroids

Forty-three lightcurves of 21 asteroids obtained in Arizona between 1968 and 1978 are presented with a brief discussion of each. Included are four asteroids not previously observed: 34 Circe, 138 Tolosa, 162 Laurentia, and 1058 Grubba. Rotation periods are at least 12 hr for Circe, either 6.42 or 12.98 hrs for Laurentia, and more than 18 hr for Grubba. Magnitudes and colors for 12 of the asteroids are given. It appears that 10 Hygiea has lightcurves which sometimes have two maxima per rotation cycle and sometimes three. A strong relation between amplitude and solar phase angle is seen for 39 Laetitia. The first direct evidence of an opposition effect for 89 Julia is given. 511 Davida is discussed in an effort to understand the pole orientation using photometric astrometry.     

Evidence for a small,high-Z,iron-like solar core

Radial and nonradial oscillation equations without and with the gravitation perturbation (with and without the Cowling approximation, CA) are solved numerically using the profile from a more accurate high-Z core (HZC) solar model. This more accurate HZC model was generated with the CRAY X-MP/48 supercomputer at the San Diego Supercomputer Center. Frequencies of oscillation in the five-min band (5MB) and frequencies with period near 160 min are presented in tables and plotted in echelle diagrams. The model was generated by integrating the stellar structure equations from the center to he surface, as done in Rouse (1964), using a maximum space step, ;x _m = 5 × 10^–4, decreasing to 10^–6 in the hydrogenionization zone just below the photosphere. Two subsets of space mesh points are used to calculate the oscillation frequencies, viz., one with a maximum space step of 5 × 10^–3, decreasing to 10^–6 with a total of 621 points (mesh 5I) and the other with a maximum space step of 2 × 10^–3, with a total of 867 points (mesh 5J).With the surface boundary condition applied at x = 1.0, the l – 1 degree nonradial frequencies with CA and the l-degree frequencies without CA are in very good agreement with the frequency spacings for observed frequencies of oscillation labeled l = 1 to 5, but with the l – 1 frequencies with CA about 10 Hz or so less than the observations and the l frequencies without CA about 10 Hz or so greater than the observations. And for the Duvall and Harvey (1983) observations labeled l = 10 and l = 20, the l = 9 and l = 19 nonradial solutions with CA agree to about 5 Hz or less with the observations. Considering from the two preceeding papers in this series that increasing the density in the outer envelope and photosphere will increase the 5MB frequencies and applying the outer boundary condition at x > 1.0 will decrease the 5MB frequencies, the net affects of such changes could move one or the other set of frequencies closer to the observations — or require a slightly different model structure to obtain accurate agreements with the values of the observed frequencies throughout the 5MB.In either case, it is concluded that the first-order, radially-symmetric structure of the model outside the HZC is close to the structure of the real Sun. This is of fundamental importance because a real gas adiabatic temperature gradient (Rouse, 1964, 1971) is used in the outer convective region without free parameters.Other aspects of agreements and differences between radial and nonradial solutions, with CA and without CA are discussed. In particular, the l = 4, 6, 8, and 9 g-mode solutions with CA indicate that the observed 160.01 min period may be a common l-mode period of oscillation. More research is proposed.